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Abstract. The yields of strange hadrons have been measured as a function of 
centrality in An + Au and in p + p collisions at ^/snn = 200 GcV in STAR. 
The system size and energy dependence are studied and compared for p+p and 
Au + Au collisions. Thermal models are fitted to the ratios of various strange 
particles to investigate the particle production and to determine the strangeness 
enhancement. The temperatures (T) and the strangeness enhancement factors 
(7s) of the systems determined from the fits are presented. 
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1. Introduction 

It was pointed out by Fermi and later by Hagedorn that particle production in 
heavy ion collisions can be described with the considerations of phase space [P. 
The question that is being discussed is whether phase space arguments can also 
describe the centrality dependence of strangeness production. The differences in the 
production of bulk matter {pT < 1 GeV) and non-bulk matter {pr > 1 GeV) are 
investigated with statistical models which are also used to estimate the equilibrium 
properties and the trends of particle yields and ratios. The quality of the statistical 
fits and the implications of the resulting fit parameters such as T, 7^ and fis are 
discussed. 



2. Statistical Models 

Assuming QGP formation in heavy ion collisions at RHIC energies, it is expected 
that the thermal nature of the partonic medium could be preserved during hadroniza- 
tion The particle yields measured in the final state then resemble a population 
in thermal equilibrium. Thermal models are used to predict the equilibrium proper- 
ties of a macroscopic system from the measured yields of the constituent particles. It 
is commonly agreed that light (u and d) quarks are more likely to reach equilibrium 
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in the medium than the strange quarks due to the relatively larger strange quark 
mass (nis ^ Tc). The strangeness saturation factor, 7^, is introduced to account for 
the amount of strangeness chemical equilibration One of the parameters of the 
statistical fits, which is the more generic term of the strangeness saturation factor 
(7s), is the phase space occupancy (7^). This term is used to regulate the sum of 
particle and anti-particle pairs produced in the medium. 

With the given ratios of particles, it is possible to deduce the temperature T, 
the volume V, the baryonic chemical potential hb, strangeness chemical potential 
Us and charge chemical potential /xq [El El El- Chemical potentials describe the 
particle anti-particle difference and require that the variables are conserved only on 
average in the whole system (Grand Canonical (GC) ensemble). The GC ensemble 
can be used to describe large systems such as central Au -I- Au collisions, while 
the Canonical ensemble is used for small system such as p + p collisions when all 
quantum numbers must be conserved exactly. 

The predicting power of statistical models for T, /is, /is, and 75 require the 
utilization of measured particle ratios. With the ratios, all degeneracy factors of 
the fireball cancel, leaving just the relative fugacities. Some examples are given in 
Equation m and Equation O 

- = K = cxp(-( )), (1) 

-=W=exp(^^^) (2) 
P T 

where the fugacity of a hadron is defined by the product of its valance quark fugac- 
ities (e.g., = A„Ad, A„ = A„A^ and Aa = A^ArfA^). 

3. Statistical Model Comparisons 

Several statistical model computational codes are available. A comparison of the 
parameters of these codes is presented in Tabled In order to make any comparisons 
of the the fits, the requirements of the models are fixed to be same. For example, 7^ 
(light-quark phase space occupancy) is fixed to 1 in SHARE and only the Grand 

Table 1. Thermal Model Comparisons [0 jT]; their ensembles and the free param- 
eters that are used in the fits. 

Models Used Ensemble Parameters 

4 Parameter Fit Grand Canonical T, ^q, fi^, and 7s 

SHARE VI. 2 Grand Canonical T, \q, As, 7,, 7s, /i/a, N, Ac and 7c 

THERMUS V2 Canonical T, B, S, Q, 7s and R 

Grand Canonical T, /is, /is, ^J.Q, fJ^c, 7s, 7c and R 
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Canonical ensemble is used in THERMUS [ 7 . Also, since the default contribution 
of the feed-down from weak decays on particle ratios is handled differently in each 
code, either all the feed-down is removed or is included in all models in the same 
amounts. 



4. Statistical Model Fits 



Predictions of a thermal model are compared to the measured particle ratios from 
p + p collisions at ^/s^ = 200 GeV in Figure ^a. The resulting fit parameters are 
Tch = 171 ± 9 MeV, 7^ = 0.53 ± 0.04 and the radius of the source, r, is 3.49 ± 0.97 
fm. The quantum numbers baryon (B) and charge (Q) are fixed to 2 while net 
strangeness (S) is fixed to 0. The short lines represent the thermal model fits from 
a canonical calculation performed by Thermus V2.0. The closed circles represent 
the particle ratios measured in p-l-p collisions. The model accurately describes both 
the stable particle and resonance ratios in p-l-p collisions. The 1 a error on the ratios 
and the differences between the model predictions and the data are included in the 
plot. The differences do not exceed 2 a. 
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Fig. 1. (a) Particle ratios in y^s^JJ = 200 GeV p + p collisions in comparison to 
canonical thermal model predictions (short lines) from Thermus V2.0 [[7]. The fit 
parameters are T^h = 171 ± 9 MeV, 7^ = 0.53 ± 0.04 and r = 3.49 ± 0.97 fm. B and 
Q are fixed to 2 while S is fixed to 0. (b) Particle ratios in the ^/snn = 200 GeV 
most central Au -I- Au collisions in comparison to a grand canonical thermal model 
prediction. The fit parameters are Tch = 168 ± 6 MeV, 7s = 0.92 ± 0.06, iib = 
(4.52±0.98) X 10-2 GeV, ns = (2.23±0.74) x 10"^ GeV, fig = (-2.05±0.77) x IQ-^ 
GeV, and r = 15 ± 10 fm. See text for details. 
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Figure njb presents the particle ratios and Thermus V2.0 predictions for the 
most central Au + Au collisions at ^^s^ = 200 GeV. The measured ratios are 
used to predict the free parameters by using a grand canonical approach. The fit 
parameters are T^h = 168 ± 6 MeV, 7^ = 0.92 ± 0.06, ^lB = (4.52 ± 0.98) x lO^^ 
GeV, us = (2.23 ± 0.74) x lO"^ QeV, fig = (-2.05 ± 0.77) x IQ-^ GeV, and 
r = 15± 10 fm. The model describes all but two particle ratios within 1 a standard 
deviation. Except for the A(1520)/A ratio, the other measured resonance ratios 
{(j>/K*, S*(1385)/A) are within a 1 a error of the model fits. Due to the very 
short lifetime (r < Tfirebaii ~ 10 fm) of most resonances, a large fraction of their 
decays occur before the thermal freeze-out. Elastic interactions of resonance decay 
products with particles in the medium alters the momenta of these particles. This 
results in a loss of the signal reconstructed for resonances (due to re-scattering). 
However, secondary interactions (regeneration) can increase the resonance yield. 
The effects of regeneration and re-scattering are not included in this thermal model 
which might explain why the experimental values differ from the statistical model 
calculations for the A(1520)/A ratio. It is expected that there should be no re- 
scattering and regeneration effects on the (j) meson due to its long life time. Since 
the model predicts the ratios for I]*(1385) and K* correctly, re-scattering appears 
to be balanced with the regeneration effects. 

In Au -I- Au collisions the free fit parameter corresponding to strangeness satu- 
ration, 7s , is higher than in p + p collisions, while T is the same within the errors. 
In Figure |3a, the energy and system size dependence of the measured K/ n ratios 
are presented An enhancement is observed in the K/tt ratios in central Au-|-Au 
collisions relative to p + p collisions implying an increase in strangeness production 
resulting in larger values of 7^ . The higher 7^ in Au + Au relative to p + p collisions 
is in agreement with the K/n ratio observations. Due to this increase of strangeness 
production from p + p to Au -I- Au collisions, particle ratios are selected such that 
/is cancels out (e.g., S*(1385)/A ratio). The small value of the chemical potential, 
fi B , can be explained by the proximity of the measured anti-particle to particle ra- 
tios (e.g., S*(1385)/I]*(1385) [ 9 ) to unity and it reflects the near zero net baryon 
number at mid-rapidity of Au + Au collisions. 

For comparisons with heavy- ion systems. Figure |21-a shows the energy and sys- 
tem dependence of the K/n ratios at mid-rapidity. Triangles correspond to data 
from heavy ion collisions and circles to p+p collisions at the given energies [ IIOL lllj . 
An enhancement in the K/tt ratios of about 50% is observed at RHIC energies in 
central Au+Au collisions relative to p+p collisions extrapolated to similar energies. 
A similar magnitude of enhancement in K~ /tt has already been observed at the 
lower energies of the AGS and SPS [d [T^. 

The energy dependence of the A and A yields at mid-rapidity from Au -I- Au 
collisions at RHIC and Pb-|-Pb collisions at SPS as a function of ^/snn is presented 
in Figure 121-b [^J. From SPS to RHIC energies, strange baryon production is 
approximately constant at mid-rapidity, whereas the A rises steeply, reaching 80% 
of the A yield at RHIC top energies. The other hyperons, S and il, follow similar 
trends. Since most of the strange baryons produced also include light up and down 
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Fig. 2. (a) Mid-rapidity K/tt ratios versus y^SNN and their dependence on heavy- 
ion and p + p collisions, (b) The energy dependence of the ratio of anti-strange to 
strange baryons from heavy-ion collisions (Pb+Pb for SPS & Au + Au for RHIC). 



quarks, strange baryon production at low energies is dominated by valence quark 
transport from the colliding system, but at RHIC it is dominated by pair production. 
This also implies that at RHIC energies, most of the incoming baryons continue 
moving with large rapidity and do not cause any increase in the baryon number in 
the mid-rapidity regions. 

5. Statistical Coalescence? 

Nuclear modification factors for strange particles in Au-I- Au collisions are presented 
in Figure ISJ-a. At high px, the ratios exhibit a suppression from binary scaling, at- 
tributed to fast moving partons losing energy as they traverse a dense medium. 
Hadron production at high pT does not follow binary scaling and is roughly con- 
sistent with participant scaling [ I15| . The clear differences between baryons and 
mesons, is believed to be due to hadron production through quark coalescence at 
intermediate pt'[I16|. For baryons and mesons, the suppression sets in at a different 
Pt- Motivated by the coalescence picture. Figure |3[b shows the Rcp ratio vs px/n 
for ^ysNN = 200 GeV, where n is the number of valence quarks. Thus pr/n rep- 
resent the Pt of a valence quark. The baryon and meson difference sets in at the 
same quark pT, in agreement with the coalescence picture. This is also observed for 
y^SNN = 62.4 GeV collisions ["TT . 

Separation of the mesons and baryons is observed when the pt of the quark 
of the meson or the baryon is in the range of 0.8-1.2 GeV which corresponds to 
a baryon pr range of 2.4-3.6 GeV and a meson pr range of 1.6-2.4 GeV. If the 



6 



S. Salur et al. 



•A 



1.8 
1.6 
1.4 
1.2 
1 

0.8 
0.6 
0.4 

.paiticipant scaling 
"■^r STAR Preliminar 



binary scaling 




? A 



5 6 7 
PT[GeV/c] 



(a) 



1.8 
1.6 
1.4 
1.2 
1 

0.8 
0.6 
0.4 
0.2 



• A 

-AK° 



binary scaling 




participant seal 



aling ? *Aa A A 



0.5 



ST^R Pre^lirtiinai^i 



1.5 2 2.5 3 3.5 

PT/n[GeV/c] 



(b) 



Fig. 3. (a) Rcp vs pT in Au + Au collisions at ^snn = 200 GeV. (b) Rcp vs 
px/n in Au + Au collisions at -^snn = 200 GeV for n = 3 for baryons and n = 2 
for mesons. Rcp is calculated from 0-5% and 40-60% central Au -I- Au collisions. 



statistical fit parameters from the ratios of the particle yields in these momentum 
ranges are compared to the ones from the integrated momentum range, it might be 
possible to determine the ranges of T and 7s for which coalescence and recombina- 
tion are applicable. The yields are extracted from the exponential distributions of 
the corrected pr spectra in the given pT ranges. The particle ratios from these sub 
Pt ranges are given in Table |21 

Table 2. The particle ratios for px of the quark 0.8-1.2 GeV that corresponds to 
baryonpT of (2.4-3.6) and meson of (1.6-2.4) in Au+Au collisions at ^/snn = 200 



Ratio 


STAR data 


Ratio 


STAR data 




1.01 ±0.02 


K-/K+ 


0.96 ±0.03 


p/p 


0.77 ±0.05 


p/n- 


0.137 ±0.013 


A/n- 


0.156 ±0.015 


A/A 


0.72 ±0.024 




0.029 ±0.003 




0.82 ±0.05 


n/n 


1.01 ±0.08 







A comparison of the T and the 7s fit parameters are presented in Figure 0] 
The Pt integrated particle ratios that are used in the fits for Figure 01-a (closed 
circles) can be found in Reference A 40% feed-down correction is applied to 

the proton yields determined from the A yields for the weak feed-down corrected 
fits (open circles). The effect of the weak feed-down correction can be seen as 
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an increase in the 7^ and a reduction in the T . For any conclusion made using 
the statistical model fits, it is essential to remember that proper weak feed-down 
contributions plays a major role in the final fit parameters. 



1.4- 

CO 

1.2- 
1- 
0.8- 
0.6- 
0.4- 
0.2- 

100 



STAR Preliminary 



"f" Kaneta et al ■■■§■■■ with feed-down 
-t-THERMUS -i- with feed-down 
• SHARE ■ -i- with feed-down 



120 



140 160 

(a) 



2.5r 



STAR Preliminary 



1.5- 



0.5- 



■■ Kaneta et al 
-THERMUS 
■ SHARE 



■4" for quark pj (0.8-1.2) 
for quark pj (0.8-1.2) 
--<>- for quark pj (0.8-1.2) 



180 



200 



Poo 1 



20 140 160 180 200 220 240 260 280 



(b) 



Fig. 4. (a) The 7^ and the T fit parameters of the integrated ratios from the given 
statistical models in comparison to the ones corrected for the weak decay feed-down. 
Input ratios of the fits can be find in [E]- (b) The 7s and the T fit parameters of 
the ratios from the quark pT range of 0.8-1.2 GeV. 

The particle ratios presented in Table |21 are used in the fits of Figure El-b (open 
circles) for comparison with the pt integrated ratios (solid circles). Since there is no 
K/tt or JI/tt ratios, the fits are less well constrained leading to larger errors as can 
be seen on the solid circles in Figure EJ-b as compared to those on the open circles 
in Figure ^a. T and 7^ increase for the quark pT range of 0.8-1.2 GeV relative to 
those for the whole pt range. Particles produced in this pT range possibly come 
from a hotter source. 



6. Conclusions 

The STAR experiment has collected tremendous amounts of data to study strangeness 
production in heavy-ion collisions at a variety of energies. The strange anti-particle 
to particle ratios show that baryon transport is approximately independent of sys- 
tem size at RHIC energies. Particle ratios in general can be used to investigate the 
fireball properties with the help of statistical models. We used several models to 
investigate the particle production in p -I- p and Au -f Au collisions at ^/snn = 200 
GeV. Comparisons of the various models and their corresponding fits are also dis- 
cussed. In these fits and model comparisons, the importance of the weak decay 
feed-down corrections is highlighted. T is similar for p + p and Au -I- Au collisions 
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at the same energy. While 7., oi p + p collisions at RHIC is smaller than the one in 
Au + Au, it is similar to that of Pb + Pb at SPS. 

Baryon and meson suppression sets in at same quark px- Since coalescence 
can partially explain the difference between baryons and mesons in y^sjviv = 200 
GeV collisions, the region where coalescence dominates the particle production can 
be investigated further with statistical models. Where quark coalescence seems to 
dominate (intermediate pr) T and 7s show an increase. This implies that particles 
in this range of pt come from a hotter source. 
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